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Abstract. Dynamics under which a system of Ising spins relaxes to a stationary 
state with Bolzmann-Gibbs measure and which do not fulfil the condition of 
detailed balance are irreversible and asymmetric. Wc revisit the problem of the 
determination of rates yielding such a stationary state for models with single-spin 
flip dynamics. We add some supplementary material to this study and confirm 
that Gibbsian irreversible Ising models exist for one and two-dimensional lattices 
but not for the three-dimensional cubic lattice. We also analyze asymmetric 
Gibbsian dynamics in the limit of infinite temperature. We finally revisit the case 
of a linear chain of spins under asymmetric conserved dynamics. 
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1. Introduction 

The one-dimensional Glauber-Ising model was probably the first example of a 
strongly interacting system with soluble dynamics showing how a system relaxes to 
equilibrium [1]. The preliminary question to solve was the choice of rates ensuring 
that a ferromagnetic chain of Ising spins relaxes towards equilibrium under single spin 
flip dynamics. The question is settled by requiring the rates at which spins flip to 
fulfil the condition of detailed balance with respect to the Hamiltonian defining the 
model [1]. 

Conversely, one may ask whether it is still possible for the system to reach a 

stationary state with the same Boltzmann-Gibbs measure by an appropriate choice 
of rates if one relinquishes this constraint, i.e., if only global balance is imposed. 
The dynamics now becomes irreversible and asymmetric: the fiipping spin is not 
equally influenced by its neighbours. This question was addressed some time ago by 
Kiinsch [2] , who exhibited examples of such a dynamics, in the particular case where 
it is totally asymmetric, in one and two dimensions. The problem was thoroughly 
revisited in ref. [3], with the following conclusions. 

The study made in [3] shows that Gibbsian irreversible Ising models exist for 
one and two-dimensional lattices but not for the three-dimensional cubic lattice. 
More precisely, imposing the up-down spin symmetry, the rate function yielding an 
irreversible Gibbsian stationary state for the linear chain depends on 3 arbitrary 
parameters. In two dimensions, the number of arbitrary parameters is respectively 
equal to 10, for the square lattice, and to 35 for the triangular lattice. Yet, for 
the totally asymmetric dynamics where only half of the spins have an influence on the 
flipping spin, the rate function is unique, up to a time scale, for these three geometries. 
In contrast, for the cubic lattice no such rate function is found, i.e., global balance 
enforces detailed balance (sec Table 4). 

The aim of the present work is to add some supplementary material to the same 
study, making the method used more easy to grasp and illustrating its outcomes on 
more examples. In particular we give a fuller account of the method, which relies on 
linear algebra coupled to the properties of the system under translation invariance, 
in order to make clearer its generality. Wc come back on the interplay between 
coordination and dimension. We explain, on the example of the linear chain, the 
constraints imposed by the positivity of the rates. We shall also be concerned in 
restating the statement made in [3], that no such rates do exist for the case of the 
three-dimensional cubic lattice, as recalled above. To this end we shall give some more 
details on the analysis in order to substantiate its conclusion. We shall then give a 
critical reading of a recent paper [4], where it is claimed that there exist irreversible 
Gibbsian dynamics for the cubic lattice, in contradiction with the study made in [3], 
and shall dismiss its conclusion on this issue. We shall finally consider the case of 
Gibbsian asymmetric dynamics at infinite temperature. An appendix is devoted to 
the study of conserved Gibbsian asymmetric dynamics for the linear chain. 

2. Stcirting point 

Let us consider Ising spins on a regular lattice of coordination in D dimensions, 
with periodic boundary conditions. The energy (Hamiltonian) of a configuration 
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C = {cTi, . . . , C7„, . . . , cTjv} reads 

E{C) = -jJ2'^n<Jj, (2.1) 

where n and j are nearest neighbours. 

2.1. Master equation 

The dynamics consists in flipping a spin, chosen at random, say spin n, with a 
rate w(C„|C), corresponding to the transition between configurations C and C„ = 
{cti, . . . , —an, • . • , CTjv}. At stationarity, the master equation expresses that losses are 
equal to gains, and reads 

P(C) wiCn\C) = MC\Cn)P{Cn)- (2.2) 

n n 

We want to find the rate function w(C„|C) satisfying this equation when P{C) is the 
Boltzmann-Gibbs distribution associated to the Hamiltonian (2.1), 

P{C) a e-^(^)/^. (2.3) 

After division of both sides by the weight P(C), eq. (2.2) can be rewritten as 

J2 {MCn\C) - w{C\CJe-^^/^) = 0, (2.4) 

n 

where the change in energy due to the flip reads 

AE = EiCn) - E{C) = 2J anhn, (2.5) 

and hn is the local field hn = J2j ^j^ "^^e to the z neighbours {dj}. We choose a rate 
function only depending on the local configuration {cr„; {ctj})} of the central spin C7„ 
and of its neighbours, and simpliiy the notation accordingly, 

w{Cn\C)=w{an;{(Tj}). (2.6) 
Thus, denoting the balance term by 

B{an; {aj}) = w{an; {a,}) - «;(-a„; {a.De-^-^""^", (2.7) 
where K = J/T, the balance equation (2.4) becomes 

YB{an;{c7j})=0. (2.8) 

n 

This equation can be satisfied either term by term, which gives the detailed balance 
condition on the rate function, 

w{<rn; {<Tj}) = w{-an; {a,})e-2^""''". (2.9) 

or as a whole, which is the global balance condition. 
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2.2. Representation of the rate function on a basis of spin operators 

The number of values taken by the rate function is equal to the number of local 
configurations {(t„: {(Tj})} of the central spin and of its neighbours. There are 2^+^ 
such configurations, i.e., 8 for the chain, 32 for the square lattice, and 128 for the 
two-dimensional triangular lattice or for the cubic lattice. We hereafter consider the 
simpler case where we have up-down spin symmetry: 

w((T„; {cTj}) = w(-(T„; {-cjj}). (2.10) 

The number of possible values of the rate function is therefore halved and is equal 
to the number of different environments of the central spin cr„, i.e., of configurations 
{(jj} of its neighbours. There are 2^ such configurations, labelled by the index a, i.e., 
4 for the chain, 16 for the square lattice, and 64 for the triangular lattice or for the 
cubic lattice. We denote the 2^ rates with a„ = -|-1 by Wa and the other 2^ rates, 
corresponding to cr„ = —1, by 

Wo, = W{an = +l\{<Jj} a), Wo, = w{an = -'\-;{(Jj}oc)- (2-11) 

The latter are obtained from the former by the spin symmetry relation (2.10), yielding 
Wa = W2- + I-CX, (2.12) 
(see Tables 1 and 2). For instance, for the linear chain, the rates to be determined are 

Wl = w{+;++), W2 = w{+:-\ — ), 

W3=w{+;-+), W4 = w{+; ). (2.13) 



Table 1. List of local configurations and corresponding values of the rate function 
for the one-dimensional chain. There are 4 possible rates Wa, with a-„ = +1, 
corresponding to the 4 possible configurations {dj}, labelled by a, of the two 
neighbours of the central spin, taicen in the order: left, right. The 4 remaining 
rates Wa, with cr„ = —1, are deduced from the former, due to the spin symmetry 
(see (2.12)). 



a 


(T„;{crj} 


Wa 


0"„ 






1 


+;++ 


Wl 




++ 


Wl = W4 


2 


+;+- 


W2 




+- 


W2 = W3 


3 


+;-+ 


W3 




-+ 


W3 = W2 


4 


+; — 


W4 






W4, = Wl 



When the spin symmetry is not imposed, the rate function depends on the values 
taken by the z + 1 spins cr^ and {uj}, and can be decomposed on a basis of 2^+^ 
spin operators made of 0, 1, . . . , 2; -|- 1 spins. For instance, for the linear chain, these 
operators are: {1, (t„_i, cr„, (t„+i, a-„(T„+i, ct„_ict„, cr„-i(T„+i, cr„_icr„(T„+i}. In the 
present situation where spin symmetry holds, this decomposition can be restricted to 
2^ even spin operators Of, i.e., 

w{crn; {cTj}) = ^ Ci Oi, (2.14) 

i=0 

with Oo — I- The knowledge of the 2^ coefficients q is equivalent to the knowledge 
of the 2^ rates Wa ■ The coefficient cq fixes the scale of time. 
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Table 2. List of local configurations and corresponding values of the rate function 
for the 2D square lattice. There are 16 possible rates Wa, with an = +1, 
corresponding to the 16 possible configurations {crj}, labelled by q, of the four 
neighbours of the central spin, taken in the order: cast, north, west, south. The 
16 remaining rates tUa, with it„ = —1, are deduced from the former, due to the 
spin symmetry (see (2.12)). 



a 




Wa 


o-„; {aj} 




1 


+ 


+ + ++ 


Wi 


- 


+ + ++ 


Wl = WiQ 


2 


+ 


+ + +- 


W2 




+ + +- 


W2 = Wi5 


3 


+ 


+ + -+ 


W3 




+ + -+ 


M)3 = Wii 


4 


+ 


+ + — 


W4 




+ H 


Wi Wi3 


5 


+ 


+ -++ 


W5 




+ -++ 


M>5 = W12 


6 


+ 


+ -+- 


Wq 




+ -+- 


= Wii 


7 


+ 


H h 


W7 




+ h 


Wl = wio 


8 


+ 


+ 


Ws 




+ 


Ws = Wg 


9 


+ 


- + ++ 


Wg 




- + ++ 


Wg = Ws 


10 


+ 


- + +- 


Wio 




- + +- 




11 


+ 


- + -+ 


Wii 




- + -+ 


wii = We 


12 


+ 


- + 


■W12 




- + 


Wi2 = W5 


13 


+ 


- - ++ 


Wl3 




h+ 


Wi3 = W4 


14 


+ 


- - +- 


wu 




- - +- 


Wii = W3 


15 


+ 


+ 


Wi5 




+ 


Wl5 = W2 


16 


+ 




W16 






W16 = Wl 



For instance, for the linear chain, 

w{an; {cTj}) = Co + Ci (TnCTn+l + C2 cr„_icr„ + C3 (T„_iCr„+i, (2.15) 

i.e., 

Ol=(7„C7„+i, O2 = Cn-iCTn, O3 = cr„_ia„+i. (2.16) 

For the square lattice, we use the following notations. The central spin fT„ being 

located at x„, we denote by aj^ (rcsp. aj^) the neighbouring spins located at x„ + ea 
(resp. x„ — Ba), where 60 (a = 1, 2) are the unit vectors spanning the square lattice. 
Thus aj^, cTjj, (Tjj and aj^ are the east, north, west and south spins, respectively. The 
list of even operators is given in Table 3. 

In order to determine the rate function satisfying the global balance 
condition (2.8), or the detailed balance condition (2.9), we can proceed in either of 
two ways. The first one consists in finding the constraints on the rates {w^}, from 
which constraints on the coefficients {cj} ensue. The second one consists in finding the 
constraints on the coefficients {ci}, from which constraints on the rates ensue. 
These two ways arc strictly equivalent because the rates {wa} arc linear combinations 
of the coefficients {cj}, and both are equivalent representations of the rate function 
wi(Tn;{(Tj}). 

We emphasize this equivalence as follows. Defining the indicator variables 

la = I{<^n = +1; Wj}a), la = I{<^n = "1; Wj}a), (2.17) 

we have, using the notation (2.11), 

w{cTn;{(Tj}) = '^{laWa+ laWa) . (2.18) 
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Table 3. List of the even operators made of the five spins ((t„; {o"j}). The central 
spin <Jn being located at Xn, <Tj^ (resp. crj^) are the neighbouring spins located 
at x„ + Ga (resp. x„ — ea), where (a = 1, 2) are the unit vectors spanning the 
square lattice. 



i 




1 




2 




3 


(T*i (7 a (T a (T a 


4 




5 




6 




7 




8 




9 




10 




11 




12 




13 




14 




15 


<^J2(^32 



Using the spin symmetry relation (2.12), we can rewrite the expression above as 

2^ 

W{an; {cTj}) = ^'Wa {la + h'+l-a) , (2.19) 
a=l 

where the two indicator variables in the bracket correspond to two opposite 
configurations. These indicator variables can be decomposed on the complete basis 
of 2^+^ spin operators made of 0, 1, . . . , 2; + 1 spins, however their sum only contains 
even operators O,. For instance, for the linear chain, 

J . f _ 1 + (T„ 1 + Sign(cr„_i)c, C7„_i 1 + sign((7„+2)„ (Tn+l 

la + h^+x-a - 2 2 

1 - C7„ 1 - sign(cr„_i)a (Jn-i 1 - sign(cr„+2) 

2 2 2 

= ^(1 + sign((7„+i)a + sign(c7„_i)a cr„_ia-„ 

+ Sign(cT„_lC7„+l)aC7„_lC7„+l) 

2^-1 

where we have introduced the matrix of signs (up to the constant 1/2^) 

A = (fli.a) = ^Oi(CT„ = +1; {aj}a). (2.21) 

(The matrix (0,^^) for the linear chain is given in the Appendix.) One can thus rewrite 
(2.19) as 

2^ 2^-1 

w{(Jn; {(Jj)) = X! X! (2-22) 

a=l i=0 
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Identifying (2.19) with (2.14) we obtain 

Ci = ^ai,a'Wa- (2.23) 
a=l 

The inverse relation reads 

2^-1 

In other words = 1/2^, where / is the unit matrix. 
2.3. Balance term 

Starting from (2.14), then using the identity e~°''^ = cosha — sinhacr, we can 
decompose the balance term on the basis of spin operators as 

2^-1 

B(a„; {a,}) = ^ 0^, (2.25) 

1=0 

where the cocfRcicnts _E,;({c,;}) arc linear combinations of the c,, with coefficients 
depending on temperature through hyperbolic functions of 2/1'. (See the Appendix 
for an illustration on the example of the linear chain.) 
Let us now define 

Ba=Wa- Wae^f^^^" , Ba=Wa- WaC^^'^" , (2.26) 

where AE^ is the change of energy associated to the rate Wq.. Thanks to the symmetry 
relation (2.12) we have 

Ba = B2^+l-a- (2.27) 
For instance, for the linear chain, (2.27) reads 

-Bi = — e^^wi = Wi — e^^Wi = B4, 

B2 = W2 - W2 = W3 - Ws = i?3, (2.28) 

and so on. Proceeding as for the rate function w{an; {o'j}), we can decompose the 
balance term B(f7„; {ctj }) as 

2^-1 

^K; Wj}) = X Fii{Ba}) Oi, (2.29) 
where the coefficients Fi{{Ba}) are linear combinations of the Ba, 

Fi = J2ai,aBa. (2.30) 

(See the Appendix for the example of the linear chain.) The sets and Fi provide 
two equivalent representations of the linear decomposition of the balance term on the 
basis of spin operators, the former expressed in terms of the coefficients q, the latter 
in terms of the values B^ taken by the balance term. 

3. Detailed balance 



We start with the simple case of detailed balance, B{an; {fj}) = 0, as a preparation 
for the sequel. This equation is satisfied by imposing Ei = Fi = for all i. 
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3.1. Constraints on the rates or on the coefficients 

The condition of detailed balance on the rate function (2.9) implies 2^ relations 
{Ba = 0}, or equivalcntly, 2^ relations {Fi = 0}. However, thanks to the symmetry 
relation (2.27) these relations are redundant and only half of them remain. We thus 
get 2^~^ relations between pairs of rates: 

w^=w^e-^''^% (a =1,..., 2^-1), (3.1) 

In return, using the spin operator representation (2.14) in these relations, the 
coefficients Ci are found to obey 2^~^ linear constraints. 

One can also proceed in reverse order, determining first the constraints on the 
coefficients Cj, then deducing those for the rates from the former. Expressing that 
-B((T„; {cj}) vanishes identically, and using (2.25) yields an homogeneous system of 2^ 
linear equations {E^ = 0}, which are not all independent. The rank of this system is 
necessarily equal to the rank of the system {Fi = 0}, i.e., to the number of relations 
between pairs of rates mentioned above, namely 2^~^. 

3.2. Examples 

We illustrate the previous considerations by the following examples. 
For the linear chain, the relations (3.1) are Bi = B2 = 0, i.e., 

w(+; ++) = e-'^^w{-; ++), 

w{+;+-)=w{-;+-). (3.2) 

The constraints on the coefficients are either deduced from (3.2) or obtained from the 
solution of the 4 equations {Ei = 0} (see Appendix A): 

ci = C2, 2ci + 7(co + C3) = 0, (3.3) 

where 

7 = tanh2ii'. (3.4) 

The space of independent rates, or independent coefficients, has dimension 2. We thus 
find the most general rate function obeying detailed balance 

w(o-„; {cTj}) = ^(1 + ^(T„_i(7„+i - |(1 + 5)cr„((T„_i +cr„+i)), (3.5) 

where a and 5 arc the free parameters, recovering a result due to Glauber, written 
here with his notations [1]. 

On the square lattice, there are 8 constraints Bi = B2 = ■ ■ ■ = Bs = 0, on the 16 
rates Wa- 

wi = e^^^ wi, 

Wa = Wa, (a = 2,3,5), 
Wa = Wa. (a = 4, 6, 7), 

W8=e^-^«)8. (3.6) 

We do not write down the corresponding 8 relations between the Cj because we will 
not use them in the sequel. 
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3.3. Symmetric rates 

Let us consider the simpler case where the rates only depend on the variation of 
energy (2.5). 

For the linear chain, the general form (3.5) automatically verifies this requirement. 
The so-called Glauber rate is obtained by fixing the parameter ^ = 0. It is the only 
such rate yielding linear equations for the temporal evolution of the observables. For 
instance, the Metropolis rate 

«;(a„;{a,}) = min(l,c-^^/^) 

^ 1 ~ TTT — {'^n-lO'n+l+ 0-nicrn-l+ O-n+l))] (3.7) 



2(1 + 7) V 2 + 7' 

does not share this property. 

For the square lattice, the requirement that the rates only depend on the variation 
of energy implies that the four neighbours of the central spin are equivalent, yielding 
11 additional constraints 

C4 = C3 = C2 = Ci, Cg = C8 = C7 = C6, 

Cl5 = Ci4 = Ci3 = C12 = Cii = Cio- (3.8) 

The remaining independent coefficients are cq, Ci, C5, ce, cio- The system of 16 
equations {Ei = 0} for these 5 coefficients only gives two constraints 

7 

Cl + -(C5 - Co) - C6 = 0, 

-^co - 2^^ci - -^C5 - Cio = 0. (3.9) 

07 D 

From this general solution one can extract some simpler expressions for the rate 
function. For instance, imposing C5 = Cio = yields 

a ( 7(2 + 7^) 
w{crn\ {<Jj}) = 2 ( ^ ~ 4(1^ ^^2) ^"(^Ji + '^h + o-ji + (^h) 

7^ \ 
+ ^^^ ^^2) ^"(^ii^J2 0'3i + '^h'^h'^n + (^h'^h'^ji + '^h'^h'^h) J(3-10) 

which is the Glauber rate, usually written as 

w{(^n; Wj}) = - (1 - C7„ tanh K{aj, + aj^ + aj, + aj^)) . (3.11) 



2 

Another simple form is obtained by setting ci = C5 = 0: 

a / _7 
2 V 2 



wicTn, Wj}) =0(1- o f^n(0-ji + (Tj^ + CT^ + (^n) 



+ ^{'^31^32 + ^^32^31 + <^JlO'i2 + ^32'^3l + ^32^32 + O-ilt^Jl)^ ■ (3-12) 

Let us finally note that the general form of a rate function satisfying detailed 
balance can be written as 

w'K; K}) = QCKDe-^"""", (3.13) 

where Q{{(yj}) is a linear combination with arbitrary coefficients of the operators O, 

not containing the central spin. This can be seen by multiplying both sides of (2.9) 
by e^'^"''" and observing that the product w(cr„; {c7j})e^'^"''" is even in (t„. 
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For instance, for the linear chain, 

<3(K}) =ao + a3 03, (3.14) 

where ao and are arbitrary, and simply related to cq and Ci of (3.3): oq = 

Co + ci tanh/C, 03 = ~(co + Ci cothJC). 
For the square lattice, 

15 

<3({CTj}) = ao + asOs + ^ a* (3.15) 

i=10 

depends on 8 arbitrary parameters. If furthermore we ask that the rate function only 
depends on the variation of energy, then one should take 

15 

Qiif^j}) = ao + a^05 + aio'^Oi, (3.16) 

i=10 

with arbitrary coefficients 00,05,010. These are linearly related to the parameters 
co,ci,C5 of (3.9). 



4. Global balance 



We now want to satisfy (2.8), not term by term but as a whole. As above we can 
solve the problem in cither of two equivalent ways: by first finding the constraints on 
the rates {Wa}, from which those on the coefficients {cj} ensue; or by first finding the 
constraints on the coefficients, from which those on the rates ensue. Both ways have 
to be implemented by formal computations, which are of equal algorithmic difficulty. 
We start with the second one, of easier presentation. 



Table 4. Number of constraints (or rank) for dynamics on regular lattices. First 
column: linear chain (2 = 2), square lattice (2 = 4), triangular lattice (2 = 6), 
cubic lattice (2 = 6). The last line is for the hexagonal lattice (see text). Second 
and third columns: number of equations and rank of the system of equations for 
detailed balance (db). Fourth and fifth columns: same for global balance (gb). 
Lasl cohiiiDi: nunibt^i' (if frcv ])araiiit*t crs Ibr .i2;l()l)al balaiK*(\ 



Lattice 


2'- 


Rank (db) 


M 


Rank (gb) 


Free parameters (gb) 


linear 


4 


2 


3 


1 


3 


square 


16 


8 


12 


6 


10 


triangular 


64 


32 


49 


29 


35 


cubic 


64 


32 


55 


32 


32 


hexagonal 


8+8 


8 


12 


8 


8 



4-1- Constraints on the coefficients 
Remind that 

B{an\{cJj])=Y,^i(^i^ (4-1) 
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Table 5. Number of constraints (or rank) and number of free parameters (gb) 
for dynamics on regular lattices at infinite temperature. Same examples as in 
Table 4. 



Lattice 


Rank (db) 


Rank (gb) 


Free parameters (gb) 


linear 


2 


1 


3 


square 


8 


6 


10 


triangular 


32 


26 


38 


cubic 


32 


29 


35 


hexagonal 


8 


5 


11 



where the Ei are linear combinations of the c,. The detailed balance condition is just 
{Ei = 0} (see section 3). The sum in (2.8) can be rewritten as 

Y,B{an;{(Tj})=NY,EiO'i, (4.2) 

n i 

defining the spatial averages of the spin operators as 

0;=^E0.- (4.3) 

n 

Taking into account the identities between the Oi due to translation invariance (see 
the examples below), the balance equation (2.8) finally reads 

^^07 = 0, (4.4) 
j 

where the Oj are a subset of the Oi, and with Eq = Eq- The size of this subset, i.e., 
the number M of terms in this sum, is equal to the difference between 2^ and the 
number of identities due to translation invariance. The M equations {Ej = 0} on the 
a are no all independent a priori. The rank of this system of equations is given in 
Table 4 for the various examples that we now present. 

4-2. Examples 

For the linear chain, we have, with the notation (2.16), 

Oi=02, (4.5) 
hence (4.4) reads 

Eo + {El + E2)Ol + EsO^ = 0. (4.6) 

The Af = 3 equations {Ej = 0} yield only one condition: Eq = (see Appendix A). 
In other words, the rank of this system of 3 equations is equal to 1 (see Table 4). We 
find: 

C1+C2+ 7(co + C3) = 0, (4.7) 

which generalizes the result (3.3) found in the detailed balance case. Hence, changing 
the notation for the coeSicients cq, ci, C2, the most general rate function satisfying the 
condition of global balance reads 

wipn, {ctj}) = -(1 - (7(1 + (5) + e)cr„c7„+i + e<T„_iC7„ 

+ 5f7„_icr„+i), (4.8) 
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which depends on the 3 arbitrary parameters a, e, S. The corresponding dynamics is 
asymmetric and irreversible. Setting e = —7(1 + S){l—p), we can alternatively write: 

+ 5f7„_iC7„+i). (4.9) 

The general Glauber form (3.5) is recovered by setting p = 1/2. At the other end, the 
case of totally asymmetric dynamics where the central spin is only influenced by one 
of its neighbours leads, once a choice of neighbour is done, to a unique expression up 
to the time scale fixed by the coefficient a. For instance, if (t„ is only influenced by 
its left neighbour, setting 6 = and p = 1 (e = —7), we obtain 

o. 

w{crn;{(Tj}) = . (4.10) 

Fixing the scale of time by the choice a = 2 cosh 2K, we obtain the exponential form 

«'K;K})=e-2^'^"-'^". (4.11) 

We shall comment further, in section 7, on the range of allowed parameters in (4.8) 
or (4.9). 

On the square lattice we have 4 identities due to translation invariance, which 
only involves two-spin operators, 

Oe = Og, Or = Os, 0^=0^3, Wi =Ou. (4.12) 
with the notations of Tabic 3. The resulting system of M = 12 (16 — 4) linear equations 
{Ej = 0}, has rank 6, yielding 6 equations of constraint on the Ci. In other words, 
6 coefflcients are expressed as linear combinations of the other 10 coefflcients, which 
remain arbitrary. Totally asymmetric cases are obtained by asking the rate function 
to depend only on two or three of the neighbouring spins instead of four. For example 
keeping the east (aj^) and north {(Jj^) spins only, and cancelling the coefficients of the 
operators containing the two other spins aj^ and aj^, we obtain a unique rate function, 
up to a scale of time, 

w{(Jn; Wj}) = f (1 - 7o-n(CTji + o-js) + 7^ o-jiC^js) • (4-13) 

Fixing this timescale by the choice a = 2cosh^2iir, allows to write (4.13) into the 
exponential form| 

w{crn; Wj}) = e-2-ff'^(^^i+<^«). (4.14) 

For the 3D cubic lattice, amongst the 2^ = 64 operators Oj, 18 are related, two 

by two, by translation invariance. They all belong to the group of (2) = 21 two- 
spin operators. The three two-spin operators of this group not related by translation 
invariance are 

for a = 1,2,3, with notations analogous to those of Table 3. So we have M = 55 
(64 — 9) linear equations {Ej = 0} in the {ci} to solve. The rank of this system of M 
equations is equal to 32, as for the case of detailed balance. The constraints are indeed 
found to be the same as when detailed balance holds (see also section 4.3 below). On 
the cubic lattice there is no irreversible Gibbsian dynamics. 



X This form, as well as (4.11), appear in [2] without the factor 2. The corresponding stationary states 
have their temperature halved. 
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It is striking to compare the former case of the 3D cubic lattice to the case of 
the 2D triangular lattice, for which the coordination is the same [z = 6). For the 
triangular lattice, there are 15 identities due to translational invariance satisfied by 
the spatial averages of a subset of the 64 operators O,. These 15 identities correspond 
to the translations of ct„(Tj^, for a = 1,2,3 (9 relations), to the translations of 
ajjCTjj , ajjCTjj ; (Tj3(Tj2 (3 relations), and to the translations of the four-spin operators, 
e.g., o-„crj3(Tji(jj2 (3 relations). Thus M = 49 (64— 15), and the rank of this system of 
equations is found to be equal to 29. In this case, there do exist irreversible Gibbsian 
dynamics. 

The totally asymmetric case involving the three spins Uj^, aj^, aj^ in the three 
unit directions is, again, determined uniquely, up to a time scale as 

w{(Tn; {cTj}) = 2 ~ 70-„(a-ji + £7j, + <Tj3) 

+7^(0-jl<^j2 + (^32(^33 + (^h'^js) - J^(^n(Tn(Tj^(Tj^) . (4.16) 

This rate function can also be written in exponential form as 

wi^n, Wj}) = e-2-^^'>(^.i+^.2+<^.3), (4.17) 
with the choice a = 2 cosh^ 2K. 



4-3. Constraints on the rates 

One can, of course, deduce the constraints on the rates from the above. Alternatively 
we can obtain these constraints directly by following the exact parallel of section 4.1. 
Remind that 

2'-l 

B{an;{(Tj})=J2 PiOi, (4.18) 

where the Fi are linear combinations of the Ba- The detailed balance condition is just 
{Fj = 0} (see section 3). The sum in (2.8) can be rewritten as 

Y,B{an;{aj})=NY,FiOl. (4.19) 

n i 

Taking into account the identities between the Oi due to translation invariance, the 
balance equation (2.8) finally reads 

J2FjO] = 0, (4.20) 

3 

where the Oj are a subset of the Oi, and with Fq = Fq. The M equations {Fj ~ 0} 
on the Ba are equivalent to the equations {Ej = 0} on the Cj. However they yield 
constraints on the rates, instead of constraints on the coefficients. 



Remark 

One can recover (4.19) following a slightly different path, as follows§. For a given fixed 
configuration C, the sum in (2.8) can be rewritten as 

[NaBc + N^B^) , (4.21) 

a=l 

§ This Vcuriant of the method was first introduced in [5]. 
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where Na (resp. No) is the number of occurrences in C of the local configuration 
{central spin up (resp. down) with the z neighbours in configuration a}, 

7V„ = 5^/K = +1;{<7,}„), 7V„ = 5^7K = -l;{ajU. (4.22) 

n n 

Using the spin symmetry and regrouping terms, the sum can be rewritten as 

Y,iNa+N2.+i-^)B^. (4.23) 

The numbers Na can be decomposed on the basis of spatial averages of all the spin 
operators. However, only the even operators remain in the sum Na + N2^+i-a 
(see (2.20)). Finahy the sum (4.23) yields (4.19). 



4-4- Examples 

We illustrate the method on the cases considered above, in one to three dimensions. 
For the linear chain, (4.20) reads 

Fo + {Fi + F2)Oi + F3O3 = 0. (4.24) 

The 3 equations {Fj = 0} yield only one condition: Fq = 0, hence the detailed balance 
condition Bi = wi — e~'^^Wi = 0, i.e., 

w{+;++)=e-^^w{-;++), (4.25) 

which is equivalent to (4.7). The rank of this system of 3 equations is equal to 1, as 
already found above (see Appendix A). 

On the square lattice, taking into account the identities (4.12) due to translation 
invariance, the resulting system of M = 12 linear equations in the B^, {Fj = 0}, has 
rank 6, yielding the following 6 equations of constraint: 

wi -e~^^wi =0, 

W6-lB6 = 0, 

W2 — e~'^^ tJ]2 + ?«5 — c~'^^ = 0, 
e^^ W3-W3- {ws - e'^^ ws) = 0, 

2 

W2 - e"^^ W2 - iw3 - Ws) + {wj - wr) = 0, 

1 + e*-^ 
2 

W2 - e~^^ W2 + - W3 - (w4 - Wi) = 0. (4.26) 

The space of independent parameters has dimension 10, in agreement with what was 
found above. 

For the 2D triangular lattice, one recovers the results presented above: the 15 
identities satisfied by the spatial averages of a subset of the 2^ = 64 operators Oj yield 
a system of M = 49 equations, the rank of which is equal to 29. 

For the 3D cubic lattice, we have M = 55 (64 — 9) linear equations in the B^, 
{a = 1, . . . , 64), to solve. The constraints found are the 32 detailed balance conditions 
Ba = 0, in agreement with what was found above. 
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Figure 1. Two types of environments for the hexagonal lattice. 



4-5. Special case of the hexagonal lattice 

On the hexagonal lattice one has to distinguish two kinds of spins, named respectively 
a and t, corresponding to two types of environments (figure 1). Spin cr„, located at 
x„, is surrounded by Tj^, Tj^ and Tj^ with notations analogous to those of Table 3, 
where ei, 62 and 63 are the unit vectors spanning the hexagonal lattice. Spin t„ = Tj-^, 
located at y„, distant from x„ by one unit, is surrounded by ak^ (which is (t„), 0-^2 
and CT/tj in the respective directions — ei, —62 and —63. To these two spins correspond 
two rates •w'^{a; {t}) and w'^{t; {a}), skipping the indices. The balance equation (4.4) 
now reads 

^ + + + ^-0j) =0, (4.27) 

j 

where the operators O" are the three o'nTj^, the three Tj^Tj^ and c^nTj-^Tj^Tj^. The 
operators O'^ are defined analogously. In (4.27) the identities due to translation 
invariance have been taken into account: 

CnTja = Tnak2, ^nJj^ — TnU k^- (4.28) 

Moreover a\7rjl and Tncrkx represent the same operator and + = counts 
for one equation only. The system of 12 resulting equations yield 8 constraints, 
identical to the detailed balance constraints. Thus, on the hexagonal lattice, there 
is no irreversible Gibbsian dynamics, i.e., global balance enforces detailed balance, as 
long as temperature is finite (see below). 

4-6. Infinite temperature 

When temperature is infinite, the balance term takes the simpler form 

-B(cr„; {(Tj}) = w{(Jn\ {(Tj}) - w(-cr„; {ctj}), (4.29) 

which only involves the operators containing the central spin. Likewise, the only 
identities between operators due to translation invariance to be considered are those 
involving the central spin. 



Rates for irreversible Gibbsian Ising models 



16 



For instance, for the linear chain, B{(Jn', {<7j}) only involves the two operators 0\ 
and 02- The identity Oi = O2 still holds, thus the balance equation (ci + C2)Oi = 
yields the constraint c\ + C2 = 0, which is the limit of (4.7) for 7 = 0. The rate 
function reads 

w{(rn\ {(J 3}) = Co + C3O3 + Cl(Oi - O2). (4.30) 

The sum of the first two terms represent the infinite-temperature reversible rate 
function Q{{(Tj}) (see (3.13)). If furthermore the dynamics is totally asymmetric, 
then one should impose ci = C3 = 0, i.e., w((T„;{(7j}) = cq, in agreement with the 
infinite-temperature limit of (4.10). In this limit the; dynamics is reversible. 

For the square lattice, i?(cr„;{aj}) only involves the 8 operators Oi to O4, and 
Oq to Og. The identities to be considered are (see (4.12)) 

0^ = 0^, 0~T = 0~i. (4.31) 

The balance equation thus imposes the 6 constraints 

Cl = C2 = C3 = C4 = 0, C6 + Cg = 0, C7 + Cg = 0. (4.32) 

This number of constraints is the same as at finite temperature. The rate function 
reads 

15 

w{cTn] {Cfj}) = Co + C5O5 + ^ CiOi 

i=10 

+ C6(06-Og) + C7(07-08). (4.33) 

As above, the first line can be identified with the infinite-temperature reversible rate 

function Q{{crj}) defined in (3.13). The second line is a linear combination of the 
operators cr„((Tj^ — CTj^) (a = 1, 2). If the dynamics is totally asymmetric, for instance 
keeping only the east and north spins, then one should impose the vanishing of the 
coefficients corresponding to operators containing the west or south spins, i.e., C5, Cg, C7 
as well as cn to C15. Thus 

w(cr„; {(Jj}) = Co -I- cioOio, (4.34) 

which corresponds to a reversible dynamics. This should be contrasted with the 
infinite-temperature limit of (4.13) which yields w(cr„;{(Tj}) = Co. There is no 
continuity of the finite-temperature result in this situation, when T ^ 00. 

For the triangular lattice, the balance term involves 32 operators containing the 
central spin. Only 6 symmetry relations remain, thus finally there are 26 constraints 
to satisfy. In other words there are more arbitrary parameters in the definition 
of the rate function satisfying global balance at infinite temperature than at finite 
temperature. As for the square lattice, the rate is equal to the sum of Q{{<Jj}) 
(32 free parameters), of a linear combination of operator differences crn{aj^ — aj^) 
(a = 1,2,3) (3 free parameters), and of a linear combination of operator differences 
of the type o"„ (criO"20"3 ^ o"30"iO"2) (3 free parameters). If the dynamics is totally 
asymmetric, for instance keeping only the spins in the direction 61,62,63, and taking 
into account the identities due to translation invariance, the resulting rate depends 
on 4 arbitrary coefficients, namely cq and the 3 coefficients eorresponcling to the 3 
operators aj^ aj^ , aj^ , aj^ aj^ . Again this dynamics is reversible but is different from 
the limit obtained from the finite-temperature result (4.16) at T — >■ 00, which yields 

w(0-„;{crj}) = Cq. 

For the cubic lattice, again the balance term involves 32 operators containing the 
central spin. Only 3 symmetry relations remain, thus there are finally 29 constraints 
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to satisfy. Detailed balance is no longer enforced by global balance at infinite 
temperature. The rate has the same form as for the triangular lattice, except that 
the four-spin operators do not enter the expression. The case of totally asymmetric 
dynamics is identical to that found for the triangular lattice. 

For the hexagonal lattice, as for the case of the cubic lattice, detailed balance is 
no longer enforced by global balance at infinite temperature. 

Table 5 summarizes the results for the various examples that we considered. 



5. Special forms of the rate function 
5.1. Totally asymmetric dynamics 

For the totally asymmetric dynamics satisfying global balance encountered so far, the 
rate function could always be written in exponential form as 

i/;K;{a,})=c-2^-"'^^, (5.1) 

where hn = + h~ is decomposed into two components related by inversion through 
the central spin see (4.11), (4.14) and (4.17). For instance on the square lattice, 

K = + f^n^ Ki = o-ji + ^n- 

Reciprocally, assume that the rate function has the form (5.1) [3]. Then the 
balance term reads 

B(a„; {c7,}) = e-2^-"''" - e-^^^"''" . (5.2) 

Expanding the two exponential terms in the right side, the analysis confirms the fact 
that on ID and 2D lattices such rates satisfy the balance equation (2.8), but shows 
that this is not the case for the 3D cubic lattice, or more generally for lattices of 
coordination ^; > 8 [3]. 



5.2. A restricted representation of the rate function 

In preparation of the discussion of section 6, we now want to investigate whether there 
are representations of the rate fimction that generalize both the form (5.1) and the 
form encoimtcrcd when detailed balance holds, namely (3.13). Let us consider, along 
the lines of Ref. [4] , the following a priori representation of the rate function, instead 
of the general expression (2.14): 

wiTn, KI) = Q(K})e-"''^", (5.3) 

where Q{{(Tj}) is a linear combination of the operators Oj not involving the central 
spin, and where 

D 

Hn = Y, (^«^Ja + Aa(TjJ , (5.4) 
a=l 

and Aa + Aa = 2K. We remind that cr^^ is the spin located at and CTj^ is the 

spin located at x„ — Ca. Setting 

Aa=K + La, Aa = K-La, (5.5) 

we can rewrite the rate function as 

w{a„;{c7,}) = Q({a^.})c-^-"''"e-^"2:.i.K„--.,J. (5.6) 

If the La vanish, one recovers the detailed balance form (3.13), made of the two first 
factors in the right side of this equation. At the other end, if La = ±K, and if Q{{<7j}) 
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is reduced to a constant, then one recovers (5.1). However thus far this form is only 

an a priori representation of the rate function, on which one should now impose the 
global balance condition. In so doing, we find the following: 

ID The form (5.3) is faithful and equivalent to the previous result (4.8). The 
constraint of global balance is automatically encoded in this form, where the 
prefactor Q{{<Jj}) = aQ + a^ O3, with arbitrary coefficients ao and 03. 

2D Let us take the example of the square lattice. The constraints of global balance 
fixes 4 linear relations on the {aj} defining Q{{(7j}) (see (3.15)). Thus the 
form (5.3) provides examples of rate functions satisfying global balance [4]. The 
reciprocal is not true, i.e., any rate function satisfying global balance is not of the 
form (5.3). For instance (5.3) implies the following relation amongst the rates: 

^=e«^^, (5.7) 

which docs not hold in general. 

One could also have argued differently by noting that the form (5.3), which 
depends on 10 arbitrary parameters a priori, actually depends on 6 parameters 
once the 4 constraints on the {aj} are imposed, namely the four remaining 
arbitrary {aj} and the two parameters Li and L2, while the general rate function 
obeying global balance depends on 10 arbitrary parameters (see section 4.2). This 
form therefore only represents a subset of the most general rate functions obeying 
global balance. 

3D Since we already stated that there is no Gibbsian irreversible dynamics for the 
cubic lattice, the form, (5.3) cannot represent a rate function satisfying global 
balance. Since the constraint of global balance enforces detailed balance, it must 
necessarily suppress the dependence of w((T„;{(Tj}) in the parameters La. In 
other words, Q{{o^j}) wist vanish for any configuration of the neighbours {cTj} 
such that (jj^ 7^ CTj^ for a = 1, 2, 3. For the configurations such that the equality 
holds, i.e., aj^ = Uj^ for all a = 1,2,3, detailed balance automatically holds, as 
it should. This is illustrated by an explicit example in section 6. 



5.3. Interpolating schemes 

For the linear chain, (4.9) gives an interpolation between the totally asymmetric cases 
and the symmetric one. 

In similar fashion there are forms of the rate function in 2D, which interpolate 
between the totally asymmetric expression (4.13) and one of the forms valid under 
detailed balance. For instance, for the square lattice, 

W{(^n; {cTj}) = 2 ~ 7PCr„((Tji + CTjJ - 7 (1 - p) an{cTj^ + (Tj^) 

+ 7^ p (Tji (Tj2 + 7^ (1 - P) <7ji CTj2 ) . (5.8) 

The totally asymmetric forms are obtained for p = 1 or 0. For p = 1/2, the rate 

function satisfies detailed balance, without being of the form (3.16). This expression 
is convenient when investigating the physical consequences of irreversibility for the 
two-dimensional Ising model with asymmetric dynamics [6]. 
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6. Review of Ref. [4] and the question of the rates for the cubic lattice 

Reference [4] is chiefly concerned with the computation of the entropy production 
rate for irreversible Ising models with Gibbsian stationary states. The rate function 
considered in this reference is of the form (5.3), which is well adapted to the 
computation of the entropy production rate. As said above, this choice is not 
restrictive in ID, but does not account for the most general rate function satisfying 
global balance for the square lattice. The expression obtained in [4] for the entropy 
production rate of the linear chain generalizes a result of [7]; the computation of 
the entropy production rate of the square lattice is done for the totally asymmetric 
case (4.14). 

We now turn to the treatment of the 3D case given in [4]. Ref. [4] claims that 
there exist irreversible Ising models on the cubic lattice, with Gibbsian stationary 
measure with respect to the Hamiltonian (2.1), contradicting the results found in 
Ref. [3] and recalled in the present work. This statement is actually untrue and relies 
on an incomplete analysis, as we now demonstrate. 

The rate chosen in [4] for the 3D case is of the form (5.3) with the particular 
choice Li = L2 = L3 = L, and where Q{{cTj}) is a linear combinations of operators, 
not involving the central spin and satisfying some additional symmetry requirements, 
with coefficients named 60, 61, . . . , 65 (see eq. (67) in [4]). It is an easy task to solve the 
problem of determining these unknown coefficients if one imposes the global balance 
condition, using the general methods described in section 4. We find the following 
constraints on the coefficients bo,bi, . . . ,b5, 

b5 = bo, 64 = 62 = —'-^-^' ^3 = ^1, (6.1) 

independently of the value of the parameter L\\. One can then check that, taking 
into account the constraints (6.1), the rate function thus obtained either satisfies the 
detailed balance condition, when (jj-^ + aj^ + Oj^ = Oj^ + Gj.^ + (Tj., , or vanishes when 
this condition does not hold, because the prefactor Q{\(yjY) vanishes itself. This is 
equivalent to saying that, for this choice of rate function, the condition of global 
balance enforces the condition of detailed balance. 

These results are in agreement with the general statement, made in [3] and in 
section 4, that there are no Gibbsian irreversible models for the cubic lattice, and 
therefore no entropy production for dynamics of the form (5.3) in this case. 

7. Positivity of the rates 

The rate functions found by the method above must satisfy the additional constraint 

of positivity for the various possible configurations. 

We illustrate the issue on the ID case, using (4.8). The allowed region in the 
plane of the two parameters (5, e) yielding positive rates is the triangle depicted 
in figure 1. The sides of the triangle correspond to the vanishing of one of the 
rates w^- The segment joining the two points (—1,-1) and (1,-7) corresponds 
to Wi = w(+;H — ) = 0. The segment joining the two points (—1,1) and (1,-7) 

II We also checked that the equations of constraint on 601 ^li ■ ■ • i ''5 written in [4] lead to the 
result (6.1). In [4] the determination of these equations of constraint is done by identification of 
two forms of the balance term, the first one deduced from the choice made for the rate the other 
one is a linear combination of a subset of the operators O, chosen according to some symmetry 
requirements. 
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Figure 2. Constraints due to positivity for the linear chain: allowed values of 
the parameters yielding positive rates and satisfying global balance are inside the 
triangle. 



corresponds to = w(+; — h) — 0. Finally the vertical segment at 6 — — 1 
corresponds to ^1=^4 = 0, i.e., w{+; ++) = w{+; ) = 0. 

All rates on the line joining (—1,0) to (1,-7) satisfy detailed balance. For 
example, the point marked M corresponds to the Metropolis rate (3.7). All rates 
with S — lead to linear equations for the temporal evolution of the observables [8] . 
The point G, located at (0, —7/2) and corresponding to the Glauber rate, is the only 
point where both detailed balance and linearity hold. 

The two ends of the green segment are the totally asymmetric points (0, —7) and 
(0,0), corresponding respectively to values of the interpolating parameter p = and 
p = 1 in (4.9). The range of allowed values with S — goes beyond this segment. 
It is comprised between the two extreme points (0, — (1 + 7)/2) and (0,(1 — 7)/2) 
which correspond respectively to the values p = (7 — l)/(27), which is negative, and 
P = il + 1)7(27), which is larger than 1. For those points and more generally for 
the range of values depicted in red in figure 1 the magnetization of the linear chain 
exhibits an oscillating relaxation [9]. For instance, for the point (0, — (1 +7)/2), the 
rate function (4.8) becomes 



w{crn;{(yj}) = 
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8. Conclusion 

The present work is a completion of [3]. One of the questions raised and solved in 
this reference concerned the possible existence of irreversible single-spin flip dynamics 
with Gibbsian stationary states for ferromagnetic Ising systems. The motivation was 
twofold. 

On the one hand, a natural question raised by the examples given in the past by 
Kiinsch [2] for totally asymmetric dynamics in one and two dimensions, is to what 
extent are these examples unique, and can they be extended to higher dimensions 
than two. The result of [3], completed here, is that, as long as the dynamics is 
not totally asymmetric, the space of parameters defining the rate function allowing 
irreversible Gibbsian Ising models is large (see Table 4). However, imposing total 
asymmetry of the dynamics yields a unique solution, up to a time scale, for the 
examples considered (linear chain, square and triangular lattices). The answer to 
the second part of the question is presumably negative. Indeed, firstly, there is no 
such Gibbsian irreversible dynamics for the cubic lattice; secondly, one can argue that 
there are neither totally asymmetric Gibbsian dynamics for lattices of coordination 
2 > 8 [3] . A novel outcome of the present work is that the situation can be different 
at infinite temperature (see Table 5). 

On the other hand, the models thus defined are interesting laboratories for the 
study of the physical consequences of irreversibility, in particular of the properties 
of the resulting nonequilibrium stationary state. For instance, for the linear chain, 
though the stationary measure is Boltzmann-Gibbs, the dynamical properties of the 
relaxing system are changed [8]. Irreversibility also implies a non- vanishing entropy 
production rate in the stationary state which can be exactly computed for irreversible 
Gibbsian models since the stationary measure is known [4, 7]. 

The method used in [3] and in the present work for the solution of the question 
raised above relies on linear algebra and properties of the system under translations. 
This implies solving the system of linear equations of constraint on the rates by a 
formal computation. It would be desirable to answer the same question by other 
means which would in some sense generalize the argument recalled above for lattices 
with coordination number ^; > 8. 
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Appendix A. Linectr chain 

Hereafter we give some details on the case of the linear chain. 
We use the identity 



(7 = tanh2/r) to obtain the decomposition (2.25) of the balance term B(cr„;{(Tj}) 
on the basis of operators Oi. The coefficients Ei thus obtained are, up to a global 
constant equal to 1/(1 — 7^), 




1 



^(l-70i)(l-702) 



(A.1) 



1 — 7' 



Eo = -7(ci + C2 + 7(co + C3)), 
El = 2ci - 7^(01 - C2) + 7(co + C3) 
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E2 = 2c2 + 7^(ci - C2) + 7(co + C3), 

= Eo, (A.2) 

from which the relation £^0/7 + (-^'i + ■£'2)/2 = is seen to hold. Hence the rank of 
the system {Ei = 0} is equal to 2. 

The matrix A = (aj^a) with < i < 3, and 1 < a < 4, which relates the 
coefficients c, to the rates Wa, reads 

^1111 

(A.3) 

Hence, using (2.30), we have, for the coefficients Fj defined in (2.29), 
Fo = {Bi + B2 + B3 + S4)/4, 
Fi = (Bi -B2+B3- B4)/4, 
F2 = (Si +B2-B3- Bi)/A, 

Fs = (Bi -B2-BS + B4)/i. (A.4) 
Noting that B2 + -B3 = 0, we obtain, up to a global constant equal to 1/4, 

Fo = (l-c4^)(,«i-e-4^Wi), 

Fi = (1 + c4^)(,i;i - e-^'^wi) - 2{w2 - W2), 

F2 = {1 + e^'^){wi - e-4^w)i) + 2(w;2 - ^^2), 

F3 = Fo, (A.5) 

thus we have the relation Fq/j + {F-^ + F2)/2 = 0. Hence the rank of the system 
{Fi = 0} is equal to 2. 

The two writings (A.2) and (A.5) can be identified by using (2.15), i.e., 

w(cr„; {cTj}) = Co+Ci UnCTn+l + C2 tT„_ia„ + C3 (7„_iCr„+i. (A. 6) 

Appendix B. Asymmetric Gibbsian conserved dynamics for the linear 
chain 

In this appendix we use the methods of the present paper to determine the rate 
function when the dynamics is conserved and satisfies global balance. So doing we 
recover the results of a previous work, which were established by a variant of the 
present method, and written differently [5]. 

Appendix B.l. Basic facts 

The dynamics of the chain consists in flipping a bond chosen at random, say bond 

(a„,(T„+i), if the two spins o-„ and (t„+i arc antialigned: either H — flips into — h, or 
— h flips into H — . The change in energy is equal to 

AF = 2 J (cr„_i(T„ + cr„+iCr„+2). (B.l) 

This is done with a rate w(ct„, (t„+i; {o-j}), where {aj} is a notation for the two 
neighbours (T„_i and an+2- The number of values taken by the rate function is 
therefore equal to 8. We denote the 4 rates with (ct„ = +l,(7„+i = —1) by Wa 
and the other 4 rates, corresponding to (cr„ = — l,cr„+i = +1), by Wa- 

Wa = w{+-]{aj}a), Wa=w{ h; {CT^jc), (B.2) 
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Table Bl. List of local configurations and corresponding values of the rate 
function for the one-dimensional chain with conserved dynamics. There arc 
4 possible rates Wa, with [an = +l,cr„+i = —1), corresponding to the 4 
possible configurations {o"j }, labelled by a, of the two neighbours of the flipping 
bond, taken in the order: left, right. The 4 remaining rates Wa correspond to 

(cTn = -l,<Tn+l = +1). 



a 


Cn, cr„+i; {cTj} 


Wa 




Wa 


1 


+-;++ 


Wi 


-+;++ 


Wl 


2 


+-;+- 


W2 


-+;+- 


W2 


3 


+-;-+ 


W3 


-+;-+ 


W3 


4 


+-; — 


W4 


-+; — 


W4 



Table B2. List of operators made of the 4 spins . . . , <Jn+2- 



i 


Oi 


1 




2 


Cn-lC„Cr„+i 


3 




4 


O'n-lO'n+lO'n+2 


5 


O'nO'n+lO'n+2 


6 


<^n-\<^n 


7 




8 


O'n-lCrn-l-2 


9 




10 




11 


Cn+lCr„_|_2 


12 




13 




14 




15 


fn+2 



(see Table Bl). 

Let us introduce the basis of 16 spin operators Oi, . . . , O15, made of the 4 spins 
(T„_i, . . . , a„+2, with Oo = 1 (see Table B2). We define the indicator variables 

la = ; {cTjla) 

la = I{ h; {cTjja) 

where Ja denotes the indicator variable of the event {{aj} in configuration a}. We 
thus have, using the notation (B.2), 

4 4 

w{+-\{<Jj}) = ^JaWa, w{- + ;{aj}) = ^JaWa- (B.4) 
a=l a=l 

Alternatively, wc have 

w(H — ; {cTj}) = Co + csOs + C12O12 + C15O15, 

«'(-+; i'^]}) = do + dgOs + di20i2 + rfisOis- (B.5) 



1 + cr„ 1 - CT. 



n+l 



2 2 
1 - cr„ 1 + (7„-|-i 

2 2 



Jcf) 



(B.3) 
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u>(c7„,c7„+i; {cTj}) = — 7r^w{+-;{aj}) 



Finally we can write the rate function as 

1 + (7n 1 - 

2 2 

1 - (7„ 1 + an +1 

2 2 

The balance term reads 



-B(o'„,cr„+i;{(Tj}) = w(cr„,(T„+i;{(Tj}) 

- w;(-<7„, {ajDe"^^/^, (B.7) 

with 

e"^''^^ = 13:^(1 -^'^6)(l-70ii), (B.8) 

(7 = tanh 2K). Using (B.5), (B.6) and (B.8), wc obtain a decomposition of the balance 
term (B.7) on the basis of operators of Table B2: 

15 

B{an,(Tr,+i;{aj}) = Eo + J2EiOi. (B.9) 
Appendix B.2. Symmetries 

We analyze the constraints induced on the rate function by the two following 
symmetries. 

Symmetry under P, the spatial left-right parity 
The constraint induced by this symmetry reads 

w{-\ — ;c7„_i,(7„+2) =w{ — h;u„+2,CT„_i), (B.IO) 

which imposes 

do = Co, ds = cg, di2 = ci5, di5 = C12, (B-H) 

or 

Wi = Wi, 'W2 = W3, W^=W2, W4 = W4. (B-12) 

Symmetry under CP 

This symmetry is the product of C and P, where the charge conjugation C changes 
the spins into their opposites. In other words, the rates are the same for + going to 
the right or for — going to the left, with the environment of the latter conjugated to 
the environment of the former. The constraints induced by this symmetry read 

w{-\ — ; a-n-i, crn+2) = w{-\ — ; -an+2, -a-„_i), 

W{ h; CTn-l, (Tn+2) = W{ h; -(Tn+2, -(Tn-l)- (B.13) 

This fixes 

C12 + ci5 = 0, di2 + di5 = 0, (B.14) 

or 



Wi = W4, Wi = W4. 



(B.15) 
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Appendix B.3. Detailed balance 



The detailed balance condition imposes -B(cr„, cr„+i; {cr^}) = 0, i.e., Ei = for all i. 

We thus obtain the 4 constraints 
1 

l~^2 



do = ^(co - 7(7C8 - ci2 + C15)), 



ds = ~ T(')'^0 + "^12 - C15)), 

rfl2 = ^ ("^12 - 7(7C15 - Co + Cs)), 

rfl5 = J3^(C15 -7(7Cl2 + co - Cg)), (B.16) 
which express the equalites 

We now restrict the rate function furthermore by symmetry requirements. 

Symmetry under P, the spatial left-right parity 

The additional constraint induced by this symmetry is 

7(co - cg) + C12 - ci5 = 0, (B.18) 
which expresses the equality 'W2 = e~'^^'W3. The resulting rate function can be read 
off from (B.5) and (B.6): 

; Wj}) = Co + C8(T„_i(T„+2 + Ci2(Tn-l + (7(co " Cg) + Ci2)o-„+2, 

W{ h; {(Jj}) = Co + C8(T„-lCT„+2 + Ci2(T„+2 + (7(co - Cg) + Ci2)cr„_i. 

(B.19) 

It depends on 3 arbitrary coefficients. 
Symmetry under CP 

The constraints on the rate function are the three first lines of (B.16) with C12+C15 = 
and di5 is fixed equal to —di2. Again the resulting rate function depends on 3 free 
coefficients. 

Appendix B.4. Global balance 

We now turn to the global balance condition. Translation invariance imposes 
0^ = 0~5, 0^ = 0^ = 0^, 0^ = 0^, 
0^2 ^O^z^ 0^= O^. (B.20) 

Solving the system of equations Ej = (see (4.4)) yields 2 constraints: 



ds. = Y3^(c8 - 7(7Co + C12 - C15)), 



7^0 - dx2 + di5 = ^-^(-ci2 + ci5 - 7(co - (2 - 7^)c8)), (B.21) 

which express the relations between rates 

W2 — W2e~^^ + — w^c^^^ = 

wi-wi+W4,-Wi + W2{1 + e"''^) - m;2(1 + e''^) = 0. (B.22) 
The number of free coefficients is equal to 6. 
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Symmetry under P, the spatial left-right parity 

This again imposes 7(00 — c%) + C12 — C15 = 0, and the resulting rate function is the 
same as for the case of detailed balance with P symmetry (see (B.19)). 



Symmetry under CP 

This imposes the two additional constraints (B.14) on (B.21), or (B.15) on (B.22). 
The resulting rate function depends on 4 free parameters. 



Appendix B.5. Totally asymmetric dynamics 

For instance, only the flipping of H — into — h is allowed. Hence w{ — h; {ctj}) = 0, or 
iBoc = 0. We thus set the two left sides of (B.21) to zero, since do = rfg = <^i2 = c^i5 = 0, 
from which it results that 

cs = 0, 7C0 + C12 - ci5 = 0, (B.23) 

or equivalently 

wi + W4 — W2 — W3 =0, W2 = wse"^^. (B.24) 

The resulting rate function reads 

w(H — ; {cTj}) = Co + ci2(T„-i + (7C0 + ci2)cr„+2, (B.25) 

which depends on 2 free coefficients. Imposing the CP symmetry fixes C12 = — 7C0/2. 
The solution found is therefore unique, up to the global time scale cq: 

wi+-,W3}) = Co (1 - |(a„_i - (Tn+2)) . (B.26) 



Partially asymmetric dynamics 

A partial asymmetry with uniform bias [10] translates into the condition 

w(-+;g„_i,cr„+2) ^ 1 - V 27) 

where <V < 1. This condition yields 

w{-\ — ; {cTj}) = Co + Ci2Cr„_i + (7C0 + Ci2)cr„+2, 
1-V 

w{ — h; {cTj}) = _^ ^ ^ (co + (7C0 + ci2)(T„-i + ci2cr„+2), (B.28) 

which depend on 3 parameters: cq,c\2 and V. Eqs. (B.24) still hold. The limiting 
case y = is included in the solution (B.19) of the fully symmetric case. The totally 
asymmetric limit V = \ reproduces the result (B.25). Imposing the CP invariance on 
the rate function again fixes C12 = — 7C0/2. 
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